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ABSTRACT 

We explore in the multi-fluid approach streaming instabilities of the electron-ion plasma 
with relativistic and ultra-relativistic cosmic rays in the background magnetic field. Cosmic 
rays can be both electrons and protons. The drift speed of cosmic rays is directed along the 
magnetic field. In equilibrium, the return current of the background plasma is taken into 
account. One-dimensional perturbations parallel to the magnetic field are considered. The 
dispersion relations are derived for transverse and longitudinal perturbations. It is shown 
that the back-reaction of magnetized cosmic rays generates a new instability with the 
growth rate much larger than that for the Bell instability, while for unmagnetized cosmic 
rays, the growth rate is analogous to the Bell's. Some difference between two models of the 
return current in equilibrium is demonstrated. For longitudinal perturbations, an instability 
is found in the case of ultra-relativistic cosmic rays. The results obtained can be applied for 
investigations of astrophysical objects such as supernova remnant shocks, galaxy clusters, 
intracluster medium, and so on. 
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1. INTRODUCTION 

It is known for a long time that a return current arises in a plasma penetrated by 
an external beam current, which is nearly equal to the imposed beam current (Roberts & 
Bennett 1968). A theory of this phenomena for the laboratory plasma has been developed 
in a number of early papers (Roberts & Bennett 1968; Cox & Bennett 1970; Hammer 
& Rostoker 1970; Lee & Sudan 1971; Berk & Pearlstein 1976). It was shown that the 
induced plasma current depends on the spatio-temporal shape of the imposed current and 
is transferred by plasma species. For external currents of the cylindrical shape, it has been 
found that the return current lies almost entirely within the beam channel (Cox & Bennett 
1970; Hammer & Rostoker 1970; Lee & Sudan 1971). In the bounded magnetized plasma 
with the given nonstationary sheet current, the return current can change with time and 
be not equal to the external current (Berk & Pearlstein 1976). However, inclusion of the 
surface current in the perfectly conducting walls results in the full compensation of both 
currents. 

The return currents in astrophysics are considered for media where cosmic rays are 
present. It is assumed that in equilibrium the total current of cosmic rays and plasma is 
equal to zero. The models are explored in which the equilibrium current is directed along 
(e.g. Achterberg 1983; Zweibel 2003; Bell 2004) and across (Riquelme & Spitkovsky 2010; 
Nekrasov & Shadmehri 2012) the background magnetic field. In the case of currents parallel 
to the magnetic field, one considers the three-component medium, the electrons, ions, and 
cosmic rays, where each species has the own drift velocity (Achterberg 1983), as well as a 
four-component one. In the last case, one assumes that the background plasma has no drift 
velocities while cosmic rays and an additional electron component (for the proton cosmic 
rays) drift together (Zweibel 2003; Bell 2004; Zweibel & Everett 2010). 

The kinetic consideration of cosmic rays drifting along the magnetic field has been 
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provided by Achterberg (1983), Zweibel (2003), Bell (2004), and Reville et al. (2007) 
for perturbations also parallel to the magnetic field. The well-known non-resonant Bell 
instability (Bell 2004) has the large growth rate for perturbation wavelengths shorter than 
the mean Larmor radius of cosmic ray protons defined by their longitudinal momentum. 
In this case, the contribution of cosmic rays to the kinetic dispersion relation is small 
(Zweibel 2003; Bell 2004). Thus, the back-reaction of cosmic rays is absent in unstable 
short-wavelength perturbations mentioned above. In the opposite case of long-wavelength 
perturbations, the perturbed currents of cosmic rays (protons) and additional electrons 
compensate each other, if only the electric drift of particles is taken into account (Zweibel 
2003). 

However, involving the polarizational current of cosmic rays in the multi-fluid approach 
is also important in the derivation of the dispersion relation for cosmic ray streaming 
instabilities. This effect was not considered e.g. by Zweibel (2003) and Bell (2004). As we 
show here, the back-reaction of magnetized cosmic rays gives rise to streaming instabilities 
which are different from that found by Bell (2004). 

In the present paper, we investigate streaming instabilities of the electron-ion plasma 
with cosmic rays up to ultra-relativistic energies in the background magnetic field, using 
the multi-fluid approach. We assume that cosmic rays, which can be both protons and 
electrons, drift along the magnetic field. One-dimensional perturbations also parallel to the 
magnetic field are treated. In this case, transverse and longitudinal movements are split. 
For generality, we take into account the thermal energy exchange between background 
electrons and ions and electron thermal conductivity. We derive dispersion relations for the 
transverse and longitudinal perturbations. For the first case, the two models with three 
and four components described above are used and corresponding results are compared. 
Analogous consideration of these models for shocks has been provided by Amato & Blasi 
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(2009). New instabilities due to the back- reaction of cosmic rays are found. 

The paper is organized as follows. Section 2 contains the fundamental equations for 
plasma, cosmic rays, and electromagnetic fields. Equilibrium state is discussed in Section 
3. In Section 4, the transverse perturbations with magnetized and unmagnetized cosmic 
rays are explored. We investigate longitudinal perturbations in Section 5. In Section 6, we 
discuss results obtained in the preceding sections. Conclusive remarks are given Section 7. 



2. BASIC EQUATION FOR PLASMA AND COSMIC RAYS 



The fundamental equations for the plasma that we consider here are the following: 

^ + v, ■ Vv, = -^L + *L E+ -^-v, x B, (1) 



dt m j n j m j m j c 



the equation of motion, 



dt 

the continuity equation, 



dfi ■ 

' 3 +v. njVj = , (2) 



dT 

_L + Vl • VT t + ( 7 - 1) T t V • = i£ (n e , T e ) (T e - T z ) (3) 



and 



dT 1 

_i + Ve . V T e + (7-l)T e V-v e = -( 7 -l)-V-q e -^(n l ,T e )(T e -T J ) (4) 

are the temperature equations for ions and electrons. In Equations (1) and (2), the subscript 
j = i,e denotes the ions and electrons, respectively. Notations in Equations (l)-(4) are the 
following: qj and rrij are the charge and mass of species j = i,e, is the hydrodynamic 
velocity, rij is the number density, pj = UjTj is the thermal pressure, Tj is the temperature, 
tf e (n e ,T e ) (iA (rij, T e )) is the frequency of the thermal energy exchange between ions 
(electrons) and electrons (ions) being vf e (n e ,T e ) = 2u ie , where v ie is the collision frequency 
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of ions with electrons (Braginskii 1965), riivf e (n e ,T e ) = n e v £ ei (rii,T e ), 7 is the ratio of the 
specific heats, E and B are the electric and magnetic fields, and c is the speed of light in 
a vacuum. We do not take into account collisions between the ions and electrons in the 
momentum equation. However, the thermal exchange should be included because it must 
be compared with the dynamical time. The value q e in Equation (4) is the electron heat 
flux (Braginskii 1965). In a weakly collisional plasma which is here considered, the electron 
Larmor radius is much smaller than the electron collisional mean free path. In this case, 
the electron thermal flux is mainly directed along the magnetic field, 

q e = - Xe b(b-V)T e , (5) 

where Xe is the electron thermal conductivity coefficient and b = B/_B is the unit vector 
along the magnetic field. We assume that the thermal flux in equilibrium is absent. 

Equations for relativistic cosmic rays we take in the form (e.g. Lontano, Bulanov & 
Koga 2002) 

+ V cr • V (RcrPcr) = h q cr E+-V cr X B , (6) 



where 



i ' cr T v c/ jrc/ / ' icr i — 1 

at n cr \ C 



In these equations, p cr = 7cr^crV cr is the momentum of a cosmic ray particle having the rest 
mass m CT and velocity v cr , q cr is the charge, p cr = % r 1 n cr T cr is the kinetic pressure, n cr is the 
number density in the laboratory frame, T cr is the adiabatic index, 7 cr = (1 — v^ r /c 2 ) ^ 2 
is the relativistic factor. The continuity equation is the same as Equation (2) for j = cr. 
Equation (8) can be used for both cold non-relativistic, T cr m cr c 2 , and hot relativistic, 
T cr ^> m cr c 2 , cosmic rays. In the first (second) case, we have T cr = 5/3 (4/3) (Lontano, 
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Bulanov & Koga 2002). The general form of the value R cr applying at any relations between 
T cr and m cr c 2 , can be found e.g. in Toepfer (1971) and Dzhavakhishvili & Tsintsadze 
(1973). 

Equations (l)-(4), (6), and (7) are solved together with Maxwell's equations 

VxE^-if (9) 



and 



„ „ 4tt. IdE 

VxB =T J+ caF' (10) 



where j = L, + j cr = £\ Qj n j v j + 3cv 



3. EQUILIBRIUM STATE 

We will consider a uniform plasma embedded in the uniform magnetic field B 
(subscript here and below denotes background parameters) directed along the z-axis. We 
assume that in equilibrium the plasma is penetrated by a uniform beam of cosmic rays 
having the uniform streaming velocity v cr0 along the z-axis. The reverse plasma current 
along this axis compensating the current of cosmic rays is provided by the streaming 
velocities of electrons, v e o, and ions, ViQ. The quasi- neutrality is satisfied due to cosmic ray 
charge neutralizaion from the background environment (Alfven 1939). Thus, we have two 
equations in equilibrium 

q e n e ov e o + qini Vio + g cr n cr0 w C ro = (11) 

and 

q e n e o + Qifiio + g cr n C ro = 0. (12) 

Such a three-component model corresponds to the one considered by Achterberg 
(1983). In papers by Zweibel (2003) and Bell (2004), the four-component model has been 
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explored in which the additional electrons (in the case of the proton cosmic rays) drift with 
the cosmic ray drift velocity. We show below that there is some difference between these 
two models. 



4. TRANSVERSE PERTURBATIONS 

We will treat one-dimensional perturbations along the background magnetic field. 
From Equations (9) and (10), it is followed that in this case the transverse and longitudinal 
perturbations are split. The transverse wave equations have the form 

The perturbed currents j lx ^ y are given by Equations (A16), (A17), (BIO), and (Bll). 
Substituting them into Equation (13), we obtain 

2 

^Ixm ~t~ t xy I1 'ix,y 



° dz 2 [dt ' ' ' 



E ix,y + £ 2 xv Ei x , y = 0, (14) 



where 



^xx ^plxx ^crxxi (15) 
&xy £plxy £crxy 

For perturbations of the form exp (ik z — iut), we find from Equation (14) the dispersion 
relation 

k 2 c 2 

^2 &xx 1 iiSxj/- (1^) 
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4.1. Magnetized Species 

We first consider Equation (16) in the case in which all species are magnetized, i.e. 

<4 » D l> ( 17 ) 

u 2 cr 3> D 2 r . 

Using Equation (17), we calculate the values e p i X x,y and e crxx ^ given by Equations (A17) 
and (Bll), respectively, and substitute them into Equation (15). Then from Equation (16), 
we derive the following dispersion relation: 

,22, X^ U pj D tj . D U lcr D lcr . ( \^ U lj D tj . W pcr Act \ , , 

where Z^^ :Cr = — iuj + ik z Vj tCr o. In Equation (18), we have neglected the contribution of the 
displacement current and small terms proportional to Dfj cr /u^ cr . According to Equations 
(11) and (12), the right-hand side of Equation (18) is equal to zero. Thus, we obtain 

a e (u - k z v e0 ) +ai(u- k z v i0 ) + a cr (u - k z v cr0 ) = k z c , (19) 

where otj = cu 2 j/cu 2 j and a cr = r y cr oR C roU 2 cr /oJ 2 cr . Solution of Equation (19) is given by 

UJ = ^ kz ±±-k z (Aj-A 1 A 3 ) 1/2 , (20) 

where 

Ai = a e + an + a cr , (21) 

A 2 = a e v e0 + a>iV i0 + a cr v cr0 , 

A 3 = a e v 2 + aivf + a cr v 2 crQ - (?. 

Using Equation (21), we find the expression A\ — AiA 3 

A\ - A 1 A 3 = A x c 2 - a e ai (v e0 - v i0 ) 2 - a e a cr (v cr0 - v e0 ) 2 - a^a^ (v cr0 - v i0 ) 2 . (22) 
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Equation (20) describes the streaming instability if (A 2 , — AiA 3 ) < 0. 

The number density of cosmic rays is considerably smaller than the density of the 
background plasma. Therefore, we can conclude from Equation (11) that v e ^ <C v cr o- In 
this case, Equation (22) can be written in the form 

A\ - A X A Z ~ (oij + a cr ) c 2 - aia cr v 2 crQ . (23) 

The growth rate of instability 5 = Imw found from Equation (20) in the case 
aia cr v 2 rQ ^> (ttj + a cr ) c 2 is equal to 

5 = (at °2 r)1/2 - k z v crQ . (24) 

This new instability arises due to the cosmic ray back-reaction, i.e. due to the same 
dynamics of cosmic rays as that of the plasma connected with the polarizational drift of the 
species. 

For the four-component model consisting of the background ions and electrons without 
drift velocities, proton cosmic rays, and additional electrons with the cosmic ray number 
density and drift velocity (Zweibel 2003), Equation (18) has the solution 

a cr j , k z r 2 21 i/2 



-k z v cr0 ± — — [-a>ia cr vl r0 + (a* + a cr ) c 2 ] , 



which is analogous to solution given by Equation (20) (see Equations (21) and (23)) except 
for the phase velocity. 



4.2. Unmagnetized Cosmic Rays 

In this section, we assume that cosmic rays are unmagnetized 

Dl » u 2 ccr . (25) 



- 10 - 



This condition can be satisfied for the relativistic cosmic rays for which 7 cr o-Rcro ^ 1- Then, 
we obtain 

U2pcr (26) 



IcroRcrO^ 2 ' 



^pcr^ccr 

&crxy 



lcroRcroD cr U 2 

The plasma ions and electrons stay magnetized. Substituting Equation (26) and e v \ xx ^ into 
Equation(16), we will have 



k 2 z c 2 -a e (oj - k z v e0 ) 2 - <x,i (uj - k z v i0 ) 2 = ±f3 cr (u - k z v cr0 ) , 



(27) 



where (3 cr = u 2 cr /u ccr . When obtaining the right-hand side of this equation, we have used 
Equations (11) and (12). We note that Equation (27) does not contain the contribution of 
the cosmic ray perturbed dynamics which is small in comparison with the plasma current 
produced by the electric drift velocities of ions and electrons. Solution of Equation (27) is 
given by 

u = — (a e k z v e0 + aik z v i0 =F J (28) 



1 

± — 

Oii 



±aif3 cr k z v cr o - a e aik 2 z (v e0 - v i0 ) 2 + -fi 2 r + aik 2 z c 2 



1/2 



From Equations (11) and (12), it is followed that v e o — ~ (?cr^cro/?i^io) v cr o (g« = —q e ). 
An estimation of the ratio of the second term in the squared brackets in Equation (28) to 
the first one gives the value (n cr0 /nj ) (k z v cr o/u ce ), which can be much less than unity. 

Solution of the dispersion relation for the four-component medium considered above is 
the following: 



_1 Per . ( , Per , Per ,22 , ^ Per , ,22^ /on\ 

u = T^— ±[ ±—k z v cr0 - - k z v cr0 + --j + k z c Al , (29) 
2 on \ on a,- w r/ . 4 



1 /2 

where c^i = (Bq / Aimiorrii) is the ion Alfven velocity and the sign || denotes an absolute 
value. We see some differences between Equations (28) and (29). 
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Equation (29) applied to the proton cosmic rays coincides with Equation (8) in the 
paper by Zweibel & Everett (2010), if we neglect the term proportional to v^ r0 (assuming 
that k z v cr0 <C | w C e I) an d take the lower sign (see also Zweibel (2003) and Bell (2004)). This 
coincidence is due to the absence of the dynamical contribution of unmagnetized cosmic 
rays to the dispersion relation. However, the nature of this effect is different in both cases. 
In our one-dimensional case, the transverse perturbations of cosmic rays does not contain 
the thermal pressure. However from the kinetic treatment, it is followed that cosmic 
rays must be cold along the magnetic field to use the transverse magnetohydrodynamical 
equations. Cosmic rays are unmagnetized under condition given by Equation (25). At 
the same time, the kinetic consideration of cosmic rays also shows that cosmic rays are 
unmagnetized in perturbations with wavelengths much smaller than their Larmor radius 
defined by the thermal velocity along the magnetic field (Zweibel 2003; Bell 2004). 

We note that if we set v £:i0 = in Equation (28) (or on the left-hand side of Equation 
(27)), we return to Equation (29) without the term ~ v^ r0 . 



We now consider potential perturbations along the background magnetic field. The 
wave equation is the following (see Equation (10)): 



In Appendices A and B, there are obtained the plasma, j p n z , and cosmic ray, j cr i z , perturbed 
currents (Equations (A18) and (B12), respectively). Substitution them into Equation (30) 



5. 



LONGITUDINAL PERTURBATIONS 



5.1. Dispersion Relation 




(30) 
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and Fourier transformation lead to the dispersion relation 



D 

T 



pit 



D 



te 



LuD te — L2 £ D ti 



q e mi 



+ 



UJ, 



pi 



Ll e Dti — L 2 iD- 



te" 



q e mi 
C[im e 



uj. 



+ -^+1, (31) 



where d/dt = —iu and d/dz = ik z . This equation will be investigated in the limiting cases. 



5.2. Cold Electrons and Ions 

We first consider the cold plasma species for which 

D% » (32) 
For cosmic rays, we here and below assume that the following condition is satisfied: 

DL » 4 r ; Tcr ° kl (33) 



where 



z? d r cr T cr o v cr0 

£/ cr — -ttcrQ n 

m cr c z c z 



In this case, the first term on the right-hand side of Equation (B13) is dominant. We note 
that the temperature of cosmic rays can be relativistic, i.e. T cr0 /m cr c 2 ^> 1. Then, using 
Equations (A6), (A8), (A10), (A12), and (B13) under conditions defined by Equations (32) 
and (33), we obtain Equation (31) in the form 

U) 2 U) 2 - UJ 2 

+ per _ „. (34) 



(u - k z v e0 ) (u - k z v i0 ) % r0 E cr0 (u - k z v cr0 ) 
where for simplicity we have neglected unity. 

Let us consider the case in which the back-reaction of cosmic rays plays the role. Then, 
solution of Equation (34) will be the following: 

w = k z v cr0 ( 1 + il~l /2 E~l /2< ^- j . (35) 



- 13 - 



If the cosmic ray back-reaction does not play the role, solution of Equation (34) will be 
given by 



The ratio of the growth rate defined by Equation (36) to that of Equation (35) is equal 
to ^l r0 E crQ (m cr /rrii) (n cr0 /nj )- Thus, the back-reaction can result in considerably larger 
growth rate even at 7 cr o 1 and T cr ^> m cr c 2 . 



Consideration shows that in the cases D\ e <C (T /m e ) k 2 z1 _D| ^> (To/m^k 2 and 
D\ <C (T /rrii)k 2 the frequency u is of the order of k z v cr0 as that in Equation (35). 
Equation (33) results in condition 7 cr0 1 when f cr0 — c. Thus, the temperature of the 
background plasma should be relativistic. However, this contradicts the basic equations 
where plasma is non-relativistic. Therefore, conditions for hot plasma are invalid. Taking 
into account other terms in Equation (B13) does not give an instability. 



We now discuss the growth rates and conditions of their derivation for transverse 
perturbations considered in Section 4. For magnetized cosmic rays obeying to Equation 
(17), the growth rate is given by Equation (24). Below, we assume that ions and cosmic 
rays are the protons. Let us consider first the case in which a-i ^> a cr or 




(36) 



5.3. 



Hot Electrons and Cold and Hot Ions 



6. DISCUSSION 



1 ^ IcroRcrO 



n crO 



(37) 
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Then, the condition of instability is the following: 



2 

IcroRcrO 3> (38) 



(see Equation (23)). The growth rate is equal to 

1/2 



I TcrO-n-crO 

V n i0 J 



k z VcrO- (39) 



This growth rate increases with the wave number k z . Let us find its value for k z = &B e ii, 
where 

_ 1 n cr Q v cr o 

^Bcll — 77 ^ci o - 

is the wave number at which there is the fastest growing mode for the Bell instability (Bell 
2004; Zweibel & Everett 2010). Then from Equation (39), we obtain 

S(k B ei\) ( D n cr0 \ 1/2 v cr0 

— 7 = IcroRcrO , (40) 

C'Bell V n i0 J CAi 



where 



1 n cr .Q v cr Q 

OBell — -W c f 



2 L-t 
^iO CAi 

is the maximal growth rate of the Bell instability (Bell 2004; Zweibel & Everett 2010). 
According to Equation (38), it is followed from Equation (40) that S (fceeii) 3> ^Beii- For the 
wave numbers such that 

k z > 



I ( 1 n cr0 \ 

™i0 / 



1/2 



CAi 

the growth rate described by Equation (39) under conditions given by Equations (37) and 
(38) is much larger then the maximal growth rate of the Bell instability. From Equations 
(20) and (21), we see that Reu <C 5. The wave numbers are limited from above by condition 
Uci/lcroRcro v cro ^ k 2 z ( see Equation (17)). From here, we can obtain an estimation of the 
maximum of the growth rate given by Equation (39) 

1/2 



j- ( 1 n cr0 \ 

Omax ~ W c i I 

\7cr0-r<cr0 UiQ J 
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The case a cr ^> a,; or 



Icr^RcrJ^- > 1 (41) 

n i0 



can be satisfied for ultra-relativistic cosmic rays for which 7 cr0 1 and/or R cr0 ^> 1. In the 
last case, the temperature of cosmic rays is relativistic one, T cr ^> m cr c 2 . The condition of 
instability has the form 

V 2 cr0 » 4" (42) 

The growth rate is given by 

5 = ( — J k z v cr0 . (43) 

\7crO-KcrO ^crO / 

A comparison of this growth rate with the Bell's one at k z = fceeii gives 

8 (fcgdl) _ / 1 ^iO \ V2 ^crO ^ 

^Bell \lcroRcrO n crO J c Ai 

The first factor on the right-hand side of Equation (44) is less than unity and the second 
one is larger than unity (see Equations (41) and (42), respectively). Therefore, the relation 
between 5 (/cbch) and <5b c ii depends on specific parameters. However, in the region 

3/2 



k z > \ {lcroRcro) 1/2 ( 

2 V / 



CAi' 



the growth rate given by Equation (43) is larger than the Bell growth rate. In the case 
under consideration, we have Reu = k z v cr o ^> S (see Equations (20) and (21)). Thus, we 
find the upper limit for k\ 

1,2 „ U li 1 H cr0 

K z ^ — 2 5 • 

V cr0 IcrO-tlcrO 

The maximal estimation of the growth rate given by Equation (43) is the following: 



IcroR 



crO 



For unmagnetized cosmic rays satisfying Equation (25) and magnetized background 
electrons and ions, solution of the Equation (27) is given by Equation (28). In the case 
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v cro/ c Ai ^ 1) the growth rate has the form ^Beii (see above). The frequency oj is much 
smaller than k-Q C \\v cr Q. Thus, Equation (25) takes the form 

TcroiicH) — 4* » !. ( 45 ) 

where we have inserted fceeii- We note that under Equation (45) cosmic rays which are cold 
in the longitudinal direction do not contribute to the dispersion relation, i.e. the cosmic 
ray back-reaction is absent. We obtain for hot cosmic rays p\\ cr ^> m cr u ccr /k z an analogous 
result, where p\\ cr is the average momentum along the magnetic field (Zweibel 2003; Bell 
2004). Substitution to the last condition the value fceeii gives 

n C rO V cr P\\ cr 



n i0 C Ai m cr 



where m rr is the rest mass. 



Let us discuss longitudinal perturbations. In the case of the cold background plasma 
expressed by Equation (32) and of Equation (33) for cosmic rays, solution of Equation (34) 
is given by Equation (35). Equation (32) can be written as v^ r0 3> T /m e . For cosmic rays, 
Equation (33) takes the form 

> —■ (46) 

Equation (46) can be satisfied for ultra-relativistic cosmic rays with v cr o — c and 7 cr o ^> 1. 
The growth rate 5 is the following: 

_3 /2 _i/2 / rn e n cr0 \ 
\m cr n e0 J 

This growth rate is much smaller than what is given, for example, by Equation (39). 

The streaming instabilities driven by cosmic rays may play a significant role in such 
environments as the shocks caused by supernova remnants (Koyama et al. 1995; Allen et al. 
1997; Tanimori et al. 1998; Vink & Laming 2003; Bell 2004), galaxy clusters (e.g. Brunetti 
et al. 2001; Pfrommer & Enfilin 2004), intracluster medium (e.g. Enfilin 2003; Guo & Oh 
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2008; Sharma et al. 2009; Sharma, Parrish & Quataert 2010) and so on, where weakly 
collisional plasma consists mainly of electrons and ions (protons) and high energy cosmic 
rays are present. Therefore, our model and results are applicable to these astrophysical 
objects. The main point of this investigation is finding that the back-reaction of magnetized 
cosmic rays can give rise to instabilities, of which the growth rate is much larger than that 
obtained earlier (Bell 2004). Although, the kinetic derivation of the dispersion relation in 
the last paper and by Zweibel (2003) contains the dynamics of cosmic rays, the contribution 
of the latter to the dispersion relation is absent in the hot regime. We obtain the same 
result for unmagnetized cosmic rays in the fluid approximation. 

7. CONCLUSION 

Using the multi-fluid approach, we have investigated streaming instabilities of the 
electron-ion plasma with relativistic and ultra-relativistic cosmic rays in the background 
magnetic field. Cosmic rays have been assumed to drift along the latter. The return current 
of the background plasma in equilibrium has been taken into account. One- dimensional 
perturbations parallel to the magnetic field have been considered. We have derived 
dispersion relations for the transverse and longitudinal perturbations, whose electric field 
is polarized across and along the magnetic field, respectively. We have shown that the 
back-reaction of magnetized cosmic rays in transverse perturbations can result in instability 
with the growth rate much larger than that for the Bell instability. For unmagnetized 
cosmic rays, we have obtained the growth rate, which is the same as the Bell's, but for 
the cold cosmic ray regime along the magnetic field when perturbation wavelengths are 
larger than the Larmor radius defined via the longitudinal momentum. For longitudinal 
perturbations, we have found an instability in the case of ultra-relativistic cosmic rays. The 
corresponding growth rate is less than that for transverse perturbations. 



The results obtained can be applied to investigation of astrophysical objects such as 
supernova remnant shocks, galaxy clusters, intracluster medium, and so on. 
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A. APPENDIX 

A.l. Perturbed Velocities of Ions and Electrons 

Let us put in Equation (1) = Vj + v^i, Pj = Pjo + Pji, E = E + E 1; B = B + Bi 
We assume that the medium and background velocities of species are uniform. Then for 
perturbations depending only on the z-coordinate and Vj || z, where z is the unit vector 
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along the z-axis, the linearized Equation (1) takes the form 

VTji TjoVwji T3 / a -i \ 

^tjVji — h -t jiH Vji x £$ , (AiJ 

mj rrijUjo rrijC 

where we have used that pji = n^T^ + n^Tjo (rij = n j0 + nji, Tj = T j0 + Tji) and 
introduced the notations D tj — d/dt + Vj d/dz and 

Fji = iL El +JZ- v x Bl (A2) 



From Equation (Al), we find equations for the perturbed transverse velocities Vji x , 



v 



i D tj + wj) Ujix = Wcj-Pjij/ + Aji^ix, (A3) 
(-°|? + ^cj) v jiy = -UcjFjix + D tj F jly , 

where u C j = qjB /mjC is the cyclotron frequency. The equation for the perturbed 
longitudinal velocity Vj iz is given by 

(5 - IS "* = ~k Dt ^ + D "^- (A4) 

where we have used the linearized continuity equation (2). 



A. 2. Perturbed Temperatures of Ions and Electrons 

From the linearized equations (3) and (4), we obtain equations for the perturbed 
temperatures of ions and electrons, T iiei . We will assume that the background ion and 
electron temperatures are equal to each other, T i0 = T e0 = T . In this case, the terms 
connected with the perturbation of thermal energy exchange frequency in Equations (3) 
and (4) will be absent. However for convenience of calculations to follow the symmetric 
contribution of ions and electrons, we formally retain different notations for the ion and 



electron temperatures. Then, we will have 



D t T a - Q te T el = - (7 - 1) T M ^, (A5) 
D e T eX - Vt ei Tn = - (7 - 1) Teo-7^- 1 . 

Here, the following notations are introduced: 

A = d u + n ie , a = D te + n x + n ei , (A6) 

^ie = v\ e (n e o, T e0 ) , fi e j = v e ei (riio, T e0 ) , 

where we have used Equation (5) for obtaining Q x . Solutions of Equation (A5) for T i>el are 
given by 

DT n = -D e (7 - 1) T l0 ^ - n ie ( 7 - 1) T e0 ^, (A7) 
DT el = -A (7 - 1) T e0 ^ - Q ei (7 - 1) T l0 ^, 

where 

-D = A -De - ftie£V (A8) 



A. 3. Equations for Longitudinal Velocities v il2 and v elz 

Let us substitute Equation (A7) into Equation (A4) written for the ions and electrons. 
Then, we obtain 

Luvaz + L 2i v elz = DD ti F ilz , (A9) 

LleVelz + L2eVi\ z = DD te F e \ z . 
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Here, we have introduced notations 





= DDI - 


- — D— — 

rrii dz 2 


(7- 


Lie 


= DDI - 


-*°D*~ 

m e dz 2 


•(7- 


Lii 


= -(7" 


L o 

1) Dti^lie 

rrii 


d 2 
dz 2 ' 


L2e 


= -(7" 


Lq 

1) D te Vlei 

m e 


d 2 

dz 2 



1^4- < A10 > 

I — AeA^, 



Solutions of Equation (A9) are the following: 



where 



Viu — ~jr (LieD ti Fii z — L 2 iD te Fei z ) , (All) 

f elz — \LliD te F elz — L 2e D u F i i z ) , 



L — L\iL\e — L2iL 2 e- (A12) 



A. 4. Expressions for Perturbed Velocities via Ei 

Using Equation (9), we can find the components of Fji given by Equation (A2). In the 
case under consideration, we obtain 

Fjix, v = ~^ Dt i i^di) Elx ' yi ( ' A13 ' ) 

-Tjlz — -friz- 

rrij 

Substitution of Equation (A13) into Equations (A3) and (All) gives 

d \ ~ l „ Qi ( d 



W + .%) ^ = ^ A, (|) " E ly + ^Dl (|) " A«, (A14) 

(di + .j) „,„ = a, (I) " A, + J^j (|) A s 
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and 



Vilz = 


D 


[ Li e D ti — 

V m i 


Velz = 


D 


[ LiiD te — 

V me 



L2iD te — Ei z , 



L2eD t i — E\ z . 

mi 



(A15) 



A. 5. Perturbed Plasma Current 

The components of the transverse perturbed plasma current j v i\ x , y = Qj n jo v jix, y are 
found by using Equation (A14). The expressions for 4n (d/dt)^ 1 j p n x , y can be given in the 
form 

-l 

jpllX = £plxxEl X + SplxyEly, (A16) 

-1 

jplly &plxyE\ x ~\~ £plxxEiy, 





d 






dt 




d 


4tt ^ 




dt 


3) 


is 



introduced in Equation (A16): 

e 



^"^f+^oUJ ' (A17) 

u&cjDtj ray 2 



The longitudinal perturbed plasma current j p n z = Ylj Qj n jo v ji + Ylj Qj n ji v jo is found 
by using Equations (2) and (A15) 

^(mT j ^^.H Ll ' D ''- Ll,D '& El ' (A18) 

,2 



pe I r n r n ^ l e 



+ -p; p- LuDt e — L 2e D t i -Eiz. 



B. APPENDIX 



B.l. Perturbed Velocity of Cosmic Rays 

The linearized version of Equation (6) for p cr i = p cr — p cr0 and v cr0 || z has the form 

^ 7 Per 1 Qcr 
RcroDfcrPcrl + VcroD tcr R cr \ = h m cr F cr iH V cr i X Bq, (Bl) 

n cr o c 

where D tcr — d/dt + v cr od/dz and 

^crF cr i = g cr Ei+^v cr0 x Bi. (B2) 

For the perturbed transverse velocities of cosmic rays v cr i Xty , we find from Equation (Bl) 
the following solutions: 



(-Dgr- + W ccr) V crlx = ^-0J ccr D tcr (^^j + ^-D cr D tcr ^J^j ^lxj (B3) 

l^cr + W ccrJ ^crly — —-^—^ccr^tcr I 7^ I -C'ls + L>cr L >tcr \ J 

where _D cr = ^ C roRcroD tcr . When obtaining Equation (B3), we have expressed F crl;E)2/ 
through E'lx^ by using Equation (A13) for j = cr (see Equation (B2)). 

The ^-component of Equation (Bl) is given by 

1 dp cr \ 

Dcr^crlz t> 'crO-RcrO-D 'tcrl "crl IcrO^crO^tcr-Rcrl o l~ Fcrlz) (B4) 



where <y crl = ^ r0 v cr0 v crlz /c 2 



B.2. Perturbed Temperature and Pressure of Cosmic Rays 

We now find R cr i and p cri . From Equation (8), we see that 
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The perturbation of the temperature T cr \ found from the equation T cr = jcrPcr/ncr is equal 
to 

rri rr I Perl n-crl . Icrl \ lT>a\ 

T cr i = T cr0 1 . (B6) 

\PcrO n crO IcrO J 

From Equation (7), we can find the pressure perturbations p cr \ 

r f n crl 7 cri \ 

Perl = Pernor , (B7) 

V^crO IcrOj 



where 

TJ.__i = TJ.._r ( 



9v cr l z /DO\ 
^crl = — n cr Q— — . \DQ) 



B.3. Equation for v cr \ z 

Substituting Equations (B5)-(B8) into Equation (B4) and using Equation (A13) for 
cosmic rays, we obtain 

3 ( t> ^ crT C rO ^ cro\ r->2 r> ^ crT C rO p, dv cr i z Y cr T cr o V cr \ z Q cr 

IcrO I K crO 5 IT ^tcr V crlz —^Icr0 v cr0 7r L> tcr — 7\ 9 — -^icr ^lz- 

\ m cr c z <r / m cr c z OZ JcrOMcr uz m c 



I'CV 



(B9) 



B.4. Perturbed Cosmic Ray Current 

The perturbed transverse cosmic ray currents j cr ix,y = qcrn C roVcrix,y are found by using 
Equation (B3). For the values An (d/dt)' 1 j cr ix,y, we will have 

4tT ( 7^ J jcrlx £ crxxE\ x ~\~ &crxyE\yi (BIO) 



d_\- 1 

dt 



4tT I J jcrly £ crxyE\ x ~\~ &crxxE\y. 
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Here 



u 2 pcr D cr D tcr ( d 



(D 2 cr + u 2 ccr ) \dt 
Up Cr uj ccr D tcr ( d x 

£ r.rxii 



(Bll) 



where u pcr = (47rn cr0 q 2 r / 'm cr ) 1 ^ 2 is the cosmic ray plasma frequency. 

The longitudinal perturbed cosmic ray current is equal to j cr i z = q C rn cr ov cr i z +q cr n cr iv cr o. 
Making use of the linearized continuity equation for cosmic rays and Equation (B9), we find 

4t (I) < B12) 



where 



, _ :'. I r> r cr T cr o V 2 t q \ 2 ^ crT C rO n 9 T cr T cr Q 2 . . 

^cr — IcrO K crO n 2~ ^tcr ~ ^IcrO^crO T^tcrT, TT^- l m <Jj 



9 9 I tcr I cru " ctu 9 — tcr o 9 ' 



